function [] = stsumplot(jptimes, fval, stimes, sval)

% STSUMPLOT Plot piecewise constant functions and their sum. The
%   functions are stored in two matrices containing the jumps points
%   and function values respectively. If they have different number of
%   jumps, each vector is padded with the last value to have the same
%   length.
%
% [] = stsumplot(jptimes, fval, stimes, sval)
%
% Inputs:
%        jptimes - a matrix of the jump times stored columnwise
%        fval - a matrix of the function values at the jump points
%          stored columnwise 
%        stimes - a column vector of the jump times of the sum
%        sval - a column vector of the function values at the jump
%          points of the sum 
%

% Authors: R.Gaigalas, I.Kaj
% v1.0 05-Oct-05

  clf;
  hold on;
  % use the built-in columwise loop
  stairs(jptimes, fval, 'b');
  % plot the sum
  stairs(stimes, sval, '-r');
  hold off;

